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ABSTRACT 

We present the first pointed X-ray observations of 10 candidate fossil galaxy groups 
and clusters. With these Suzaku observations, we determine global temperatures and 
bolometric X-ray luminosities of the intracluster medium (ICM) out to rsoo for six 
systems in our sample. The remaining four systems show signs of significant contam¬ 
ination from non-ICM sources. For the six objects with successfully determined rsoo 
properties, we measure global temperatures in the range 2.8 ^ Tx ^ 5-3 keV, bolomet¬ 
ric X-ray luminosities of 0.8 x 10"'^'^ ^ Tx,boi ^ 7-7 x erg s“^, and estimate masses, 
as derived from Tx, of M 500 ^ 10^^ M 0 . Fossil cluster scaling relations are constructed 
for a sample that combines our Suzaku observed fossils with fossils in the literature. 
Using measurements of global X-ray luminosity, temperature, optical luminosity, and 
velocity dispersion, scaling relations for the fossil sample are then compared with a 
control sample of non-fossil systems. We find the fits of our fossil cluster scaling re¬ 
lations are consistent with the relations for normal groups and clusters, indicating 
fossil clusters have global ICM X-ray properties similar to those of comparable mass 
non-fossil systems. 

Key words: X-rays: galaxies: clusters - galaxies: clusters: general - galaxies: groups: 
general 


1 INTRODUCTION 

Fossil galaxy systems are group and cluster mass objects 
characterized by extended, relaxed X-ray isophotes and an 
extreme magnitude gap in the bright end of the optical lu¬ 
minosity function of their member galaxies. Typically, fos- 
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sils are identified with the criteria of a halo luminosity of 
7/x,boi ^ 0.5 X 10"^^ erg and a first ranked galaxy more 
than 2 i?-band magnitudes brighter than the second bright - 
est galaxy within half the virial radius ll Jones et al.ll200^ . 
Fossil systems comprise 8-20 per cent o f groups and clus - 
ters in the same X-ray luminosity regime (I Jones et al.ll2003l . 
and thus determining the origin of the features characteriz¬ 
ing these systems is important for understanding the nature 
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and evolution of a significant fraction gala:xy groups and 
clusters. 

The features of fossil systems seem to fulfil theoreti¬ 
cal predictions that the Milky Way luminosity (L*) galax¬ 
ies in a group will merge into a central bright elliptical in 
less than a Hubble time, but the time-scale for the cool¬ 
ing and collapse of the ho t gas halo is longer l|Barneslll989l : 
IPonman fc Bertramlll99^ 1. Indeed the first fossil group dis¬ 
covered, RX J1340.6+4018 (IPonman et al.l[i994l i. appeared 
as a solitary bright elliptical located in the centre of a group¬ 
sized X-ray luminous halo. It was thought the central galaxy 
of this group was the final merger remnant of the former 
group galaxies, and hence this object was named a ‘fos¬ 
sil group’. Since then, deeper observations have found this 
system to con sist of galaxies ot her than the bright central 
galaxy (BCG; Ijones et al.ll200(Ji and as a result the magni¬ 
tude gap criterion of fossils has been established. The mo¬ 
tivation for this criterion is that over time, an increasingly 
growing difference between the two brightest galaxies will 
form as a result of the merging of the most massive galax¬ 
ies into a single bright central elliptical if no infall occurs. 
This formation scenario is well suited for group mass fos¬ 
sils where the velocity dispersion is low and the dynamical 
friction time-scale is short. 

A number of objects meeting the fossil criteria 
have als o been observed in t he cluster mass r egime 
as well j jCvDriano et all l2006l : iKhosroshahi et H] l2006l : 
Voeyodkin et al. 120101: lAguerri et al.l 1201 ll : iHarrison et al.l 


2 OI 2 II . It is possible fossil clusters may form as the result 


of two systems merging, where one group has had its bright 
galaxies merge due to dynami cal friction, and the o ther has 
comparatively fainter galaxies dHarrison et al.ll2012l l. Should 
merging occur between systems with similarly bright galax¬ 
ies, any previously existing magnitude gaps may become 
filled in. Therefore, meeting the fossil criteria may only be 
a transitory phase in the evolution of a group or cl uster 
dvon Benda,-Beckmann et al. l l2008l : [Dariush et al.ll201o| ). 

Numerical and hydrodynamic simulations indicate the 
large magnitude gaps characterizing fossil groups and clus¬ 
ters are associated with an early formation time: fossil 
systems have been found to assemble more of their total 
dynamical mass tha n non-fossil systems at every redshift 
dPariush et al.ll2007l ) , where half the dyna mical mass is as¬ 
sembled hy z > 1 dD’Onghia et al.| I2OO5I I. Evidence that 
fossils have formed and evolved in a different manner than 
normal groups and clusters should then manifest in differ¬ 
ences in their respective properties. 

The bright central galaxy which dominates the optical 
output of fossil systems has a number of unique character¬ 
istics, although whether this demonstrates a clearly distinct 
formation scenario from non-fossil BCGs is still uncertain. 
The BCGs of fossils are more massive in both the stellar 
component and in total than the cen tral ellipticals in non ¬ 
fossil systems of the s ame h alo mass (IHarrison et al.l[2012ll . 
iMendez-Abreu et al.l d2012l l find fossil BCGs are consistent 
with the Fundamental Plane of non-fossil BCGs, but show 
lower velocity dispersions and higher effective radii when 
compared to non-fossil intermediate-mass elliptical BCGs of 
the same A'g-band luminosity. These results suggest the fos¬ 
sil BCG has experienced a merger history of early gas-rich 
dissipational mergers, followed by gas-poor dissipationless 
mergers later. 


On the global scale, the scaling relations of fossil sys¬ 
tems remain a point of contentio n due to limited da t a and 
inhomogeneities between studies. IKhosroshahi et ^ ll2007l . 
hereafter KPJ07) performed a comprehensive analysis of a 
sample of group mass fossil systems and found their sample 
fell on the same Lx-Tx relation as non-fossils. However, the 
fossil groups were found to have offset Lx and Tx for a given 
optical luminosity Lopt or velocity dispersion (j„ when com¬ 
pared to normal groups, which was interpreted as an excess 
in the X-ray prope rties of fossil systems for their mass. In a 
comparable study, IProctor et al.l 1 I 2 OIII I found similar devi¬ 
ations between fossils and non-fossils. This offset, however, 
was interpreted as fossils being underluminous in the opti¬ 
cal which is supported by their large mass-to-light ratios. 
These features would not result from galaxy-galaxy merging 
in systems with normal luminosity functions, and thus this 
analysis calls into question the formation scenario commonly 
attributed to generating the characteristic l arge magnitude 
gap o f foss il systems. Later stu dies, such as IHarrison et al.l 
ll2(il2tl and iGirardi et aP (l2014l . hereafter G14), find no dif¬ 
ference in the Lx-Lopt relation of fossil systems and non - 
fossils. Even so, most recently IKhosroshahi et ^ (l2014h 
present a sample of groups, one of which qualifies as a fossil, 
that lies above the Lx-Lopt relation of non-fossil systems, 
reopening the debate on fossil system scaling relations. 

In this paper we have undertaken an X-ray study of 
10 candidate fossil systems, never previously studied with 
detailed pointed observations in the X-ray regime. Using 
Suzaku data, we present the first measurements of intra¬ 
cluster medium (ICM) temperatures, bolometric X-ray lu¬ 
minosities, and estimates of the M 500 masses of our sys¬ 
tems. This work comprises the sixth instalment of the FOs- 
sil Group Origins (FOGO) serie s. The FOGO proj ect is 
a multiwavelength study of the ISantos et al. (1200711 can¬ 
didat e fossil system catalogue. In FOGO I ( Aguerri et al.l 
I2OIIII . the FOGO project is described in detail and the 
specific goals of the collab oration are outlined. FOGO H 
dMendez-Abreu et alJ I 2 OI 2 I I presents a study of the BGG 
scaling relations of fossil systems and the implications for 
the BCG merger history. Global optical lumi nositi es of 
our FOGO sample are measured in FOGO HI llG14h and 
used to construct the global Lx-Lopt relation which re¬ 
veals no difference between the fossil and non-fossil fits. 
Deep r-band observations and an extensive spectroscopic 
database were used to redetermine the magnitude gaps of 
the FOGO sample and reclassify our fos sil candidate cata¬ 
logue in F OGO IV (IZara ttini et al.ll2014l . hereafter Z14). In 
FOGO V (IZarattini et al.ll2015ll . the correlation of the size 
of the magnitude gap and the shape of the luminosity func¬ 
tion is investigated. In this work (FOGO VI) we advance the 
characterization of the X-ray properties of fossil systems and 
constrain the global scaling relations of these objects. 

The details and observations of our Suzaku sample are 
described in Sections [5] and |31 A discussion on how non- 
ICM sources may contribute to the observed emission of our 
systems follows in Section [I] Tests to determine the contri¬ 
bution of these non-IGM sources are presented in Sections [S] 
and ini Measurements of the global ICM properties of the 
thermally dominated subset of our sample are recorded in 
Section [T] Global scaling relations and their implications are 
presented in Section[8l For our analysis, we assume a ACDM 
cosmology with a Hubble parameter Ho=70 km s”*^ Mpc“^, 
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a dark energy density parameter of Qa=0.7, and a matter 
density parameter r2M=0.3. 


2 THE SAMPLE 


Our sample of 1 0 observed ga l ajcy g roups and clusters was se¬ 
lected from the lSantos et alJ (12007, he reafter S07) catalogue 
of candidate fossil systems. The IsOtI catalogue was assem¬ 
bled by first identifying luminous r <19 m ag red galaxies in 
the lu minous red galaxy (LRG) catalogue dEisenstein et al] 
l200iri . and selecting only those galaxies associated with 
extended X-ray emission in the ROSAT All-Sky Survey 
(RASS). Sloan Digital Sky Survey (SDSS) Data Release 5 
was then used to spatially identify companion galaxies to 
these bright galaxies. Group or cluster membership was as¬ 
signed to galaxies identified within a radius of 0.5 Mpc 
from one of the bright LRGs and with a redshift consistent 
with that of the LRG. While spectroscopic redshifts were 
used when available, galaxy membership was primarily de¬ 
termined using photometric redshifts. Groups and clusters 
with more than a 2 r-band magnitude difference between 
the brightest and second brightest member galaxies within 
the fixed 0.5 Mpc system radius were then selected, and 
tho se wi th an early-typ e BC G were identified as fossils. 

\Zl4 observed the ISOTI fossil candidate list with the 
Nordic Optical Telescope, the Isaac Newton Telescope, and 
the Telescopio Nazionale Galileo to obtain deeper r-band im¬ 
ages and spectroscopic redshifts for candidate group mem¬ 
bers allowing for improved system membership. Addition¬ 
ally, the search radius for galaxy system members was ex¬ 
tended to the virial radius of th e syst em as calculated from 
the RASS X-r ay lu minosity. The lZ14l study confirms 15 tar¬ 
gets out of 34 [sot] candidates are fossil galaxy systems. Ac¬ 
cording to this characterization, our sample contains five 
confirmed fossil systems and five non-confirmed or rejected 
fossil systems (see Table [T]). 


3 OBSERVATIONS AND DATA REDUCTION 

The 10 systems in our sample were observed with the 
Suzaku X-ray telescope between 2012 May and October (Ta¬ 
ble [T|. Our analysis uses the data from Suzaku's three X- 
ray Imaging Spectrometers (XIS) sensitive to the 0.5-10 
keV band. Our single-pointing observations were taken with 
a normal clocking mode, and an editing mode of 3x3 or 
5x5 which were combined when both were available. The 
stacked XIS0+XIS1+XIS3 raw count images of the sample 
are shown in Fig. [T] 

The analysis of our study was conducted using the 
HEASOFT version 6.15 software library with the calibra¬ 
tion database CALDB XIS update version 20140520. Spec¬ 
tra were extracted using XSELECT version 2.4c and fit us¬ 
ing XSPEC version 12.8.Ig. The event files were reprocessed 
using aepipeline with the CALDB XIS update 20140203 
using the default settings with an additional criterion of 
COR>6. In our spectral analysis, emission from the ®®Fe 
calibration sources, located in the corners of each XIS de¬ 
tector, was removed. Additionally, the XISO damaged pixel 
columns caused by micro-meteorites were masked. 

A Redistribution Matrix File (RMF) was created for all 


spectral extraction regions with xisrmfgen. For each RMF, 
two Ancillary Response Files (ARFs) were created with 
xissimarfgen, one to be convolved with the background 
spectral model, and the other to b e convolved w i th th e 
source model following the method of llshisaki et al.l ll2007tl . 
Background ARFs were created out to a radius of 20 arcmin 
using a uniform emission source mode. For the source ARFs, 
an image of the stacked XIS field-of-view (FOV) was used 
to model the emission. 


4 TREATMENT OF NON-ICM EMISSION 

High fidelity measurements of the IGM temperature and lu¬ 
minosity require careful consideration of non-IGM sources 
of emission during our analysis. 


4.1 Background and foreground sources 


The standard Suzaku XIS back ground consists of a non-X- 
ray particle background (NXB; Tawa^.et_al, 2008ll. th e cos¬ 
mic X-ray background lGXB: [Fabian fc Barconslll992ll . and 
foreground Galactic emission from the Local Hot Bubble 
(LHB ) and the Milky Way Halo fMWH: lKuntz fc Snowd^ 

I2OOOII . 


The contribution of the NXB for each object was as¬ 
sessed using the night earth database within 150 days of 
the o bservation using the FTOOL xisnxbgen (iTawa et al.l 
I2OO8II . Our XISl observations were taken in a charge injec¬ 
tion mode of GI = 6 keV which increases the NXB. Accord¬ 
ingly, the nxbscie calibration file was used as input for 
XISl to counteract this. 

The contribution of the galactic foreground to a 
XIS spectrum is well described by two thermal plasma 
models: apeclhb+(wabs x apecmwh) where 2:lhb = 
2mwh = 0, Zlhb = Zmwh = 1 Zq, and feTcHB = 
0.1 keV (iKuntz fc SnowdenI I2OO0I I. The CXB was mod¬ 
elled by an absorbe d power-law: wabs x powerlawcxb 
with r = 1.412 dKushino et al.l l2002h . During spec¬ 

tral analyses, the summed background and foreground 
model: apeclhb+wabs(apecmwh+powerlawcxb) was 
convolved with the uniform emission ARF. 


4.2 Solar wind charge exchange 

The interaction of ions in the solar wind with neutral atoms 
in the heliosphere and in Earth’s atmosphe re can produce 
E < 1 keV photons in the X-ray regime llGravensI 120001 : 
iFuiimoto et aLlboOTlL To check for contamination from so¬ 
lar wind charge exchange (SWGX), proton flux light curves 
with a sampling frequency of 90 s were obtained from the 
NASA WIND-SWE database over the time span of each ob¬ 
servation. The intensity of proton flux has been found to be 
related to the strength of geocoronal SWGX contaminating 
photons, where flux levels above 4 x 10® protons cm“^ s“^ 
commonly indicate potentially signific ant contamination t o 
X-ray spe ctra from charge exch ange (lYoshino et al.ll2009ll . 
Following IFuiimoto et all (l2007l l. 2700 s were added to the 
time points in the WIND-SWE light curve to account for 
the travel time between the WIND satellite, located at the 
LI point, and Earth, where the geocoronal SWGX emission 
is produced. 
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Table 1. Summary of observations. 


Object 

Sequence number 

RA 

Dec. 

Start date 

Exposure |ks] 

TypeT 

FGS03 

807052010 

07;52;44.2 

4-45:56:57.4 

2012 Oct 28 18:39:14 

14.3 

F 

FGS04 

807053010 

08:07:30.8 

4-34:00:41.6 

2012 May 06 16:24:20 

10.1 

NG 

FGS09 

807050010 

10:43:02.6 

4-00:54:18.3 

2012 May 30 05:18:38 

9.9 

NG 

FGS14 

807055010 

11:46:47.6 

4-09:52:28.2 

2012 May 29 17:06:08 

12.4 

F 

FGS15 

807057010 

11;48;03.8 

4-56:54:25.6 

2012 May 26 17:58:41 

13.6 

NF 

FGS24 

807058010 

15:33:44.1 

4-03:36:57.5 

2012 Jul 28 08:10:10 

13.2 

NF 

FGS25 

807049010 

15:39:50.8 

4-30:43:04.0 

2012 Jul 28 18:06:02 

10.6 

NF 

FGS26 

807054010 

15:48:55.9 

4-08:50:44.4 

2012 Jul 29 02:05:54 

8.6 

F 

FGS27 

807056010 

16;14;31.1 

4-26:43:50.4 

2012 Aug 05 07:14:36 

10.6 

F 

FGS30 

807051010 

17:18:11.9 

4-56:39:56.1 

2012 May 02 11:43:31 

14.0 

F 


t The fossil status column contains the lZ14l updated fossil characterizations of the lSOTl catalogue. In the fossil 
status column, ‘F’ is a co nfirm ed fossil, ‘NF’ is a rejected fossil, and ‘NC’ is not confirmed as either a fossil 
or non-fossil according to lZ14l and remains a fossil candidate. 


Much of the FGS24 observation occurs during an el¬ 
evated period of proton flux; however, the light curve of 
FGS24 displays no significant duration flares. Furthermore, 
as a check, we have performed our spectral analysis on 
the time windows where the proton flux was less than 
4 X 10®cm“^s“'^ and found the results were consistent with 
the spectral analysis of the full baseline. We therefore con¬ 
sider the effects of SWCX to be small and have recorded the 
results of the analysis of the full observation in the main text 
and include the FGS24 light curve and shortened exposure 
time analysis in Appendix lAl 


4.3 Point source contamination 

Our Suzaku observations are the first pointed X-ray ob¬ 
servations of the objects in our sample. Consequently we 
must assess point source contamination primarily relying on 
the Suzaku data alone. Visual inspection of the XIS images 
(Fig. H]) reveal two obvious point sources in the FGS15 FOV 
which we are able to exclude in our analysis using circular re¬ 
gions of radius 2.5 arcmin. Additionally, FGS03 and FGS09 
show diffraction spikes from a strong point-like sources near 
the peak of the X-ray emission. However, the large 2 arcmin 
half-p ower diameter (HPD) o f the Suzaku X-ray Telescope 
fXRT: [Serlemitsos et al.ll2007fl inhibits the exclusion of these 
sources and the robust identification of other point sources. 

Optical and radio studies of the objects in our sample 
have found a number of active galactic nuclei (AGN) in spa¬ 
tial proximity to our galaxy systems. Especially concerning 
are the radio-loud AGN, located near the projected location 
of the BCGs, fou nd in 7 out of the 10 objects in our sample 
l|Hess et al ]mB)- To determine if these radio-loud AGN, 
and other optical and radio AGN in the FOV, are signifi¬ 
cant contributors to the source emission in the X-ray regime, 
we perform image (Section [5]) and spectral (Section |6]) anal¬ 
yses. In the 0.5-10 keV range of the XIS, the strength of 
AGN emission increases towards the harder energies of the 
spectrum. As a result, the harder photons from an AGN 
may falsely boost the measured temperature of the ICM if 
only a thermal model is used to fit the spectrum. Assessing 
AGN contribution is therefore a crucial step in determining 
the properties of the IGM. 


4.4 Implementation of RASS data 

Because most of our objects extend over the entire sin¬ 
gle Suzaku pointing, a local Suzaku background region is 
not consistently available to assess the background contam¬ 
ination in our source regions. To aid in constraining the 
LHB, MWH, and CXB, we employ RASS background spec¬ 
tra sensitive to the 0.1-2.4 keV X-ray regime. RASS spec¬ 
tra were obtained through the High Energy Astrophysics 
Science Archive Research Center (HEASARC) X-ray back¬ 
ground tool0 in an annulus of inner radius 0.5 degrees and 
outer radius 1 degree centred on each of our sources. The 
size of this annulus is sufficient to minimize contamination 
from the source itself where the largest rsoo radius found for 
an object in our sample only extends to ~20 per cent of the 
inner radius of the annular RASS background region. 


5 IMAGE ANALYSIS 

5.1 Determination of the source aperture for the 
spectral analysis region 

The region of our initial spectral analysis for each object was 
established to encircle where the emission from the source 
dominates the emission from the background, enabling the 
parameters describing the source spectrum to be determined 
in a high signal-to-noise ratio (S/N) region. We determine 
this source region using vignetting and exposure corrected 
images of the source as well as simulated images of the back¬ 
ground estimated from RASS spectra. 

For each Suzaku pointing, an exposure map was created 
with xisexpmapgen and a flat-field using xissim. The flat- 
held was simulated over the XIS 0.5-10 keV energy range at 
a monochromatic photon energy of 1 keV for a uniform sky 
out to 20 arcmin. 

An image of the NXB particle background for each 
pointing was produced with xisnxbgen over the same energy 
range. This image was estimated from night Earth observa¬ 
tions within 150 days of the Suzaku observation date. The 
NXB image was uniformly corrected by dividing the count 
rates by the exposure time. 

Emission from the CXB, LHB, and MWH 
was estimated from RASS background spectra. 

^ http;//heasarc.gsfc.nasa.gov/cgi-bin/Tools/xraybg/xraybg.pi 
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Figure 1. The Suzaku combined raw counts XIS0+XIS1+XIS3 images in the 0.5— 
10 keV band. The image is Gaussian smoothed with a = 0.42 arcmin. White circles 
demarcate the initial spectral extraction region rap,src defined to encircle the source- 
dominated region (rap,src values in Table[4||. ^^Fe calibration source events have been 
removed. 


These spectra were fit with the background model: 
apeclhb+wabs(apecmwh+powerlawcxb). Because the 
RASS background spectrum consists of only 7 data points, 
only the normalizations of the three background compo¬ 
nents were allowed to vary; the other parameters were fixed 


at the standard literature values as described in Section HD 
The ROSAT PSPC response matrix provided by the 
background tool was implemented for the fit. In calculating 
the background photon flux in the Suzaku XIS 0.5-10 keV 
energy range, the XSPEC dummyrsp command was used to 
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Table 2. General information. 


FGS“ 

Coordinates of Peak X-ray ^ 

2 = 



RA 

Dec. 


[1020cm-2] s 

03* 

07:52:46.48 

+45:56:48.40 

0.052 

5.06 

04 

08:07:29.47 

+34:01:02.95 

0.208 

4.27 B 

09 

10:43:03.33 

+00:54:33.26 

0.125 

3.88 1 

14* 

11:46:47.37 

+09:52:33.38 

0.221 

2.89 1 

15 

11:48:02.43 

+56:54:49.57 

0.105 

0.998 s 

24 

15:33:43.74 

+03:37:03.74 

0.293 

3.65 1 

25 

15:39:49.57 

+30:42:58.40 

0.097 

2.29 

26* 

15:48:56.03 

+08:50:51.27 

0.072 

3.14 

27* 

16:14:30.77 

+26:44:02.18 

0.184 

3.61 C 

30* 

17:18:11.79 

+56:39:51.33 

0.114 

2.21 1 


^ [SMS2007] ID ^ 

^ Coordinates determined from the stacked XISO+XISl+XISS s 

raw count image in the 0.5-10 keV band I 

^ Spectro scopi c redshift of the central bright galaxy in the fossil 
cluster <|S0^ | 

^ Weighted average gala ctic hydrogen colum n density in the 
direction of the target llKalberla et al.l|2005h 


FGS30 


10'^“ 


10 ® 


10 '’ 


10 ® 

10 '® 


lO'io 



radius (arcsec) 


radius (arcsec) 


Confirmed fossil system 


Figure 2. An example of the estimated source and background 
surface brightness profiles for FGS30. The bottom right-hand 
panel shows the average source and background profile for the 
three XIS detectors. 


extrapolate beyond the ROSAT PSPC sensitivity range of 
0.1-2.4 keV. 

An image of the estimated CXB+LHB+MWH emission 
was produced with xissim out to a radius of 20 arcmin from 
the coordinates of the X-ray centre of the systems. The emis¬ 
sion was modeled with the best-fitting spectral model and 
photon flux of the RASS background data. Because of the 
low count rate of CXB+LHB+MWH photons over the ex¬ 
posure time for each object, the exposure time was increased 
by a factor of 10, and corrected later, to improve the statis¬ 
tics of the surface brigh tness profile of the result ing image 
following the method of iKawaharada et al.l (120101 ). 

An image of the source could then be created from 
the images constructed during this procedure. Because the 
NXB background is not affected by vignetting, the expo¬ 
sure corrected image of the NXB was subtracted from the 
exposure corrected image of the XIS detector. The result¬ 
ing image was then vignetting corrected with the flat-held 
and the vignetting and exposure corrected image of the 
CXB+LHB+MWH was subtracted to obtain the estimated 
vignetting corrected image of source emission. 

Surface brightness profiles were created using ds9 for the 
vignetting corrected source, NXB, and CXB+LHB+MWH 
images as shown for example in Fig. [5] The coordinates of 
peak X-ray emission (Table [2|) were used as the centre of 
the surface brightness profile. The profile waa constructed 
from 20 uniformly spaced circular annuli out to the radius 
of the largest circle that could be inscribed within the XIS 
FOV from the centre coordinates. The source and combined 
background profiles were then averaged for the three XIS 
detectors and the radius at which the source and background 
emission are equal was identified. We find that within this 
radius the source contributes on average -^80 per cent of 
the total counts, with no less than a ~70 per cent source 
contribution for all objects in our sample. It is this radius, 
the source radius r^p.src, which we have used to define our 
region of initial source spectral analysis. 


5.2 Surface brightness analysis 

Radial surface brightness profiles were constructed for each 
object using stacked 0.5-10 keV XISO+XISl+XISS observed 
images. For the purpose of this profile analysis, we apply 
an additional satellite attitude correction to the event files 
used to create the images. Suzaku XIS images can contain 
up to a 1 arcmin position error as a result of a recurrent off- 
set between the XRT optical axis and the satellite attitude 
llUchivama et al.ll2008l ). With the application of a corrected 
attitude file, the XIS images can thus be sharpened. This 
correction was performed by generating corrected attitude 
files with aeattcor, and then applying these corrected at¬ 
titude files to our cleaned event files using xiscoord. The 
new corrected event files are used to produce the images 
used in our brightness profile analysis, the brightness pro¬ 
files of which are shown in Fig. [3] The number of annuli 
for each profile was determined such that each annulus had 
at minimum 225 counts, which, assuming Poissonian noise, 
requires the number of counts to be 15 times the error. 

The brightness profile of a spherically symmetric and 
isother mal ICM in hydrostatic equ i librium will follow a /?- 
model (ICavaliere fc Fusco-FemiaM 1 19761 . Il978l ) . These are 
appropriate assumptions for virialized and relaxed groups 
and clusters. Disparity between the data and the sin¬ 
gle /?-model can therefore result from processes such as 
merger asymmetries, multiple thermal components, and 
non-thermal emission, for example, as produced by an AGN. 
Our initial fit of the profiles consists of a /3-model plus a 
background constant: 

5'(r) =5'o(l + (r/rc)")-®^+^/" + fc, (1) 

where So is the central surface brightness, Tc is the core 
radius, and k is the background surface brightness. In this 
model, the /?-model component was convolved with a ra¬ 
dial model of the Suzaku XRT PSF (see Appendix [0. Fits 
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Table 3. Best-fitting parameters of the surface brightness pro¬ 
files. 


result, we compare the fit of three source models to our 
spectra: 


^-model + background constant 



cT 

*0 

rc 

P 

k^ 

xVd.o.f. (x/ 


[10-2| 

|kpc| 


[10-^1 


03^ 

299-9l37;3 


1 - 00 ^ 0.05 

89.2+j_| 

90/27 (3.3) 

04 

26.7t”, 

48tl 


6 7+°-2 
' - 0.2 

25/15 (1.6) 

09 

50.4ltj 

3Stt 

l.OOj^Q Q5 

19 9 + 0-3 

126/22 (5.7) 

14* 


90+ 1 1 

^■^-'^-0.02 

4 3+0-® 

65/30 (2.2) 

15 


letl 

0 49 +®-°^ 

26.31“:^ 

99/19 (5.2) 

24 

9 7+0.8 

38+j? 

0 40 '*'^'°^ 

Q 3 

3.9+K:| 

49/28 (1.8) 

25 

34 4+2'® 


0 45 + 0.01 

o-olU 

80/40 (2.0) 

26* 

12 7+^'® 
-0.9 

47tl 

0 37+0-01 
-0.00 

0.0+!;® 

23/17 (1.4) 

27* 

a-olKI 


U.OO_Q Q5 

9 8+°-^ 
^•®- 0.4 

37/21 (1.8) 

30* 

so.olj®;? 

1©^ 

0 39+0-01 
0-Oy Q QQ 

0 4 + 2-5 
^■^- 2.7 

60/40 (1.5) 


^ Units of counts s ^ Mpc ^ 
Confirmed fossil system 


were performed with the Sherpa Python module (iDoe et all 
1200711 . 

The returned best-fitting parameters are recorded in 
Table |3] and the convolved best-fitting model is shown in 
Fig. 131 We note that FGS03, FGS09, FGS15 have Xr > ^ 
indicating the /1-model poorly describes the observed emis¬ 
sion. For these objects, we test adding to the original model 
a point-like component consisting of a 5 function convolved 
with the PSF model. This additional point-source compo¬ 
nent does not offer an improvement in Xr compared to the 
original /3-model fits. Nevertheless, the emission from these 
three objects seems to indicate that either the ICM is not 
relaxed, or there is some significant source of non-ICM emis¬ 
sion. 

Because the annuli used are smaller than the Suzaku 
XRT PSF and, additionally, discrepancy from a /3-model 
could be attributed to multiple phenomena, we consider the 
results as merely suggestive and to be used and interpreted 
in conjunction with our spectral analysis. 


6 SPECTRAL ANALYSIS 

Our spectral analysis consists of measuring spectral proper¬ 
ties within a region of high S/N lSection l6.ll) and using these 
results to classify these objects as thermally dominated or 
AGN contaminated lSection l6.2ll . The results of this section 
will then be used to measure or estimate the global proper¬ 
ties of the IGM-dominated systems within rsoo (Section [7|. 


6.1 Spectral fitting in the source region 

In order to disentangle ICM emission from potential con¬ 
taminating point source emission, we perform our analysis 
on the source aperture region where the source emission is 
more than half of the total emission from the object. By de¬ 
termining this source aperture radius, rap,arc as described in 
Section we make no assumptions on the type of source 
emission. Extracting a spectrum from this region therefore 
improves the spectral analysis of any type of source over 
the background whether the source is dominated by ther¬ 
mal emission from the ICM or non-thermal emission from 
an AGN. 

The results of our surface brightness profile analysis in¬ 
dicate some objects in our sample may have a strong non- 
thermal point-like component to the total emission. As a 


(i) an absorbed thermal plasma model, wabsxapec, to 
model the ICM; 

(ii) an absorbed power-law, wabs x powerlaw, to model 
an AGN; 

(iii) an absorbed combined thermal and power-law model, 
wabs(apec+powerlaw), to describe contribution 
from both the ICM and an AGN; 

where the wabs absorption component accounts for galactic 
absorption in all three models. 

The background and foreground sources consist of the 
NXB, LHB, MWH, and CXB. The NXB spectrum was used 
as the background file for the extracted Tap region to be sub¬ 
tracted directly during the spectral fit. The CXB, LHB, and 
MWH were accounted for through modelling as described in 
Section \n\ 

The XIS spectra were grouped with grppha such that 
each bin had a minimum of 25 counts. The binned Suzaku 
XISO, XISl, and XIS3 spectra were fit simultaneously with 
the RASS background spectrum. The Suzaku spectra were 
fit with the source and background model while the RASS 
spectra were fit only with the background model. The RASS 
best-fitting parameters were tied to that of the Suzaku spec¬ 
tra with a scaling factor to account for the difference in the 
angular size of the spectral extraction regions. Bad channels 
were ignored for all spectra. The Suzaku XISO and XIS3 
spectra were fit over 0.7-10 keV fSection 16.1.Ill , the XISl 
spectra over 0.7-7 keV, and the RASS spectra over the range 
0.1-2.4 keV. 

In all three models, the neutral hydrogen column den¬ 
sity was assigned the w eighted average galac tic value in the 
direction of the source llKalberla et al.ir2005l l. The redshifts 
of our systems were taken to be the spectrosco pic r edshifts 
of the bright central galaxies as determined bv ISOTI . During 
the fit, the column density and redshift were always fixed. 
The metal abundance Z component of t he apec model was 
calcul ated using the abundance tables of lAnders &; Grevessd 
lll989ll . The photon index of the po werlaw model was con¬ 
strain ed to be within F = 1.5 — 2.5 jlshibashi fc Gourvoisieil 
l20ldl . Initially, all other parameters were left free to be fit. 
However, if during the fit convergence on an apec or pow¬ 
erlaw parameter within the physically reasonable limits did 
not occur or the parameter was returned with infinite error 
bars, the fit was performed again with that parameter fixed. 
In all further tables, quantities presented without error bars 
have been fixed to a reasonable value. 

The resulting best-fitting parameters are listed in Ta¬ 
ble |4] and the best-fitting models to the spectra are shown 
in Fig. 131 The background parameters resulting from each 
of the model fits were consistent with each other within la 
errors. 

6.J.1 A soft energy excess 

While the XIS is sensitive to photons with energy as low as 
0.5 keV, we have excluded the E < 0.7 keV energy channels 
from our spectral analysis. In the majority of our observa¬ 
tions, an apparent excess in counts was found in the 0.5-0.7 
keV range when compared to the fit of the apec or power- 
law models in the E > 0.7 keV range. 
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radius [Mpc] radius [Mpc] radius [Mpc] 

FGS30 



radius [Mpc] 


Figure 3. Surface brightness profiles of the stacked XIS image in the 0.5-10 
keV band. The best-fitting convolved /3-model is plotted in solid red; dashed lines 
represent the components to the model. Residuals for the /3-model are plotted as 
squares. 


Potential origins of this soft excess include a second 
thermal component in the ICM, an AGN, calibration issues, 
SWCX, or statistical fluctuations. Adding a second thermal 
model to the ICM model did not improve the fit. If an AGN 
were the origin of the excess, removing the softest energies 


should not greatly deter detecting the presence of its emis¬ 
sion in the spectra because an AGN will contribute most 
strongly to the harder energies of the spectrum. Calibration 
issues with proportional removal of flickering pixels from ob¬ 
servations of the source and the NXB may also contribute 
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to energy channels below 0.6 keV. Additionally, it is possi¬ 
ble there is some contribution from SWCX in the soft energy 
regime, although the solar wind proton flux light curves of 
most of our sample are of a low intensity indicating geocoro¬ 
nal SWCX emission is unlikely to be a significant contami¬ 
nant (see Section [4.2|l . 

Because the origin of this excess is uncertain and thus 
cannot be appropriately modelled in the spectra, and fur¬ 
thermore the excess only affects the first few low energy 
channels in the spectrum, we exclude this softest energy 
regime from our fits. This has little effect on the returned 
best-fitting parameters and in general the reduced of the 
fits improves with the exclusion of the soft excess energy 
channels. 


6.2 Comparison and interpretation of the model 
fits 

In comparing the fits of the three models, the 
APEC-fPOWERLAW model does not appear to signifi¬ 
cantly improve the characterization of the spectra over 
the individual apec and powerlaw fits. Indeed in the 
combined fit, the apec and powerlaw components are 
never simultaneously constrained. As a result, while some 
APEC-fPOWERLAW fits return with values slightly less 
than that for the less complex fits of apec or powerlaw 
only, we decide to choose the simpler model that has all 
parameters constrained. 

By the Xr values, the powerlaw model provides a bet¬ 
ter fit over the thermal apec model for FGS03, FGS09, 
FGS15, and FGS24. We consider these four objects to be 
dominated by a non-IGM source and with our current ob¬ 
servations, we cannot disentangle the IGM and non-ICM 
emission. Further discussion on these objects is provided in 
Appendix El 

For the remainder of our analysis, we focus on those 
objects in our sample with spectra that are best fit by the 
APEC model and are thus galaxy systems dominated by IGM 
emission. 


7 GLOBAL ICM TEMPERATURES AND 
LUMINOSITIES 


In order to compare the ICM temperatures and luminosities 
of our fossil systems with those of other groups and clusters, 
we calculate these properties within the fiducial radius of 
rsoo, the radius at which the average enclosed density is 
500 times the critical density of the Universe. We calculate 
rsoo, and the spectral properties within this radius, using an 
iterative procedure. 

Using the temperature calculated within some aper- 
ture. Tap, we calcul ate rsoo using the rsoo-Tx relation of 
lArnaud et al.l ll2005ll : 


rsoo = 1.104 h'jQ E(zy 


kT V 


5 keV, 


Mpc, 


( 2 ) 


where /i 7 o = Hq/{70 km s ^ Mpc and E lz) = 

H(z)/Ho = V^m{ 1 + zy + IIk(l + zy + flA (lHogelil999ll . 
This value of rsoo is used as our next radius of extraction to 
determine a new Tap, and we continue this process until con¬ 
vergence is reached between rsoo and the temperature, and 
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FGS 14 (apec) 



ifi TO'i'ji II' I'T 
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Energy (keV) 
FGS 25 (apec) 
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Energy (keV) 

FGS 15 (powerlaw) 



1 2 5 


Energy (keV) 
FGS 26 (apec) 



"i 2 ^ ^ 5 

Energy (keV) 

FGS 24 (powerlaw) 



1 2 5 


Energy (keV) 
FGS 27 (apec) 




Energy (keV) 




Figure 4. The XISO (black), XISl (red), and XIS3 (green) spectra for the source regions 
rap,src determined in Section 15.11 The best-fitting model to the observed spectra, as 
determined by the values in Table (4] is plotted in a solid line. The RASS spectra 
that were simultaneously fit with the Suzaku background model are not shown. 


thus Tsoo has been determined. This analysis is performed 
on the subset of our sample that is thermally dominated 
(Section 16.2[) . The iterative process is begun with the Tap 
determined from the apec only fit as recorded in Table [1] 

For two of our objects, FGS25 and FGS26, the final 
estimation of rsoo extends beyond the largest aperture ra¬ 
dius that can be inscribed within the XIS FOV. However, 
our estimated rsoo is very similar to the largest aperture 
size that was used to extract spectral parameters, where 
the ratio between the maximum rap and rsoo is 0.98 and 
0.84 for FGS25 and FGS26, respectively. As a result the Tap 
values for these two objects should reasonably describe the 
true global temperature within rsoo- When considering the 
luminosity, Tx, 5 oo is estimated from Tx,ap using a surface 


brightness profile model that well describes the IGM emis¬ 
sion. By integrating this surface brightness model over area, 
the conversion factor between Tx ,500 and Tx,ap is calculated 
using the relation 

T^x.soo ^ S{r)r dr 
ix,ap S{r)r dr 

where S is an azimuthally averaged radial surface brightness 
profile. For our surface brightness model, we use the /?-model 
where So, r^, and /3 have the values recorded in Table [31 
With the global temperature values listed in Table (Sj 
we estimate the masses within rsoo for our systems using 
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the Msoq-Tx relation of lArnaud et alj ll2005l '): 

Msoo = 3.84 X 10^" E{z)-^ (^) 

We find our thermally dominated objects have masses con¬ 
sistent with clusters (M 500 ^ 10^“^ Mq). 


8 SCALING RELATIONS 

We combine our newly measured global Z/x,boi ,500 and Tx 
with previously measured fossil systems properties, to con¬ 
strain the scaling relations of these objects with the goal of 
assessing if fossil systems display different scaling relations 
than those for normal groups and clusters. Our analysis of 
fossil system scaling relations is distinguished from previous 
studies through several updates including the fitting of the 
largest assembled fossil system data set, using recent X-ray 
and optical data for our control sample of normal groups and 
clusters, and a substantial effort of homogenizing both the 
fossil and non-fossil data sets. We furthermore record the 
best-fitting Lx-Lr relation, and for the first time record the 
slopes and y-intercepts of the Lx-7x, Lx-ffv, Ty^-a^ scaling 
relation fits for fossil systems. 


8.1 Sample assembly, correction, and fitting 


We have assembled data from a number of studies to inves¬ 
tigate how the global X-ray and optical properties of fossil 
systems compare to non-fossil groups and clusters. To en¬ 
sure a reliable comparison, we have made an effort to use 
quantities determined within the same fiducial radius and 
defined the same way. For our analysis we use bolometric X- 
ray luminosities Lx.boi, temperatures Tx, and optical SDSS 
r-band luminosities Lr all calculated within rsoo, and global 
velocity dispersions (t„. While we have removed known fos¬ 
sils from our sample of non-fossil groups and clusters, we 
do not have information on the magnitude gap between the 
first and second brightest galaxies in all of the systems mak¬ 
ing up our control sample. However, the large magnitude 
gap characterizing fossil systems should be fou nd in only a 
fract i on of Lx.boi ^ 5 x 10"^^ erg s“^ systems (Ijones et al.l 
l2003l : lMilosavlievic et al1l2006ll . Thus, we expect our control 
sample is contaminated by at most a few unidentified fossil 
systems. 

To assembl e our group sample, we us e the (j„ of the 
‘G-sample’ from lOs mond &: Ponmari ||2004 ). Gr o up Tx val¬ 
ues ar e pulled fromiRasmussen fc PonmanTi 2007|) , Sim et ahl 
j2009|}j_lHudsonet_^ (l201(]l l , lEckmiller et al.l ll201ll l , and 


iLovisari et al. 


iet_^ |2010|}j_^ckmillCTetalJ ll201lll . 

1 2 OI 5 II . ILovisari et al. ( 2015l l Lx , 0 . 1 - 2 . 


are transformed to L x.boi using the conversion tables of 
iBohringer et al.l ll2004ll . 

For the cluster sample, we use the r-band op¬ 

tical luminosities calculated within rs op. The corre s pond¬ 
ing v e locity dispersions are take n fromiGi rardi et al.l lll998l . 
200^ Girardi fc Mezzettil (l200ll l. |Popesso et al. 1 200lf l , and 


Zhang et al.l 1 20 fil l. Bolometric X-ray luminosities within 


rsoo and temperatur e s were sourced from IPratt et al.l ll2009l l 
and iMaughan et al. ll2012ll. and su pplemented with addi¬ 
tional Lx.b oi from Zhang_et__alJ ) 201ll l and Tx from lWu et al.l 
||1999D an d iHudson et al.l ll20ld h 

Taking our sample of fossil systems observed with 


Suzaku, we match the global X-ray properties of t he sys¬ 
tems in T able [ 5 ] with the correspon ding Lr from lGl3 and (7„ 
from IZl j . For the remainder of the lZldl confir med fossil cat¬ 
alogue, we supplement the Lx.boi from 1^. For improved 
consistency with the Lx of our cluster sample, we approx¬ 
imate X-ray luminosities that more closely resemble those 
computed using the gr owth curv e ana lysis (GGA) method 
dBohringer et al.l l2000f l from the IG14| lu minosities derived 
from RASS counts (see section 3.3 of lGldl for details). These 
corrected luminosities also show good agreement with the 
Suzaku measur ed Lx for the sample of objects shared be¬ 
tween both the lGldl study and ours. 

We add to the fossil s ample w ith the X-ray luminos i- 


ties and temperatures from IKPJOTI and Millei_e£_^ 


matched with the Lr and (t„ data from IProctor et ^ 


( 2012 ), 


( 201111 . 


The lKPJOTl Lx.boi,200 are rescaled to rsoo using their best¬ 
fitting /3-model parameters and our luminosity conversion 
Eq. [31 To ensu re consistency with our Suzaku sample, the 
rsoo of lKPJO^ is recalculated from their temperatures using 
our Eq. [ 2 ] and we use this value t o esti mate Lx.boi.soo- To 
rescale the Lr. .200 of IProctor et al.l (1201 il l to rsoo, we assume 
the light follows the mass, which is a go od approximation 
on th e global scale of groups and clusters dBahcall fc Kulien 
I 2 OI 4 II . For a NFW density profile w ith concentration pa - 
rameter c = 5, M500/M200 = 0.70 dNavarro et al.l [ToOTII . 
The assumption of c = 5 was chosen for agreement with 
the concentrati ons of normal clusters of similar tempera¬ 
ture and mass (|Pointecouteaii_etjiy[200^; |Pratd_j£^^Ammi^ 


200E; IVikhlinin et al.l 20061 : Buote et al.l 20071 : Ettori et al 


2 OIOII because the typical concentration parameter for fossil 


systems is poorly characterized. Thus, we can rescale us¬ 
ing Lopt, 500 /Lopt ,200 oc 0.70. Here, the correction relation is 
applied only to the non-BCG light. _ 

W e also implement the fossil catalogue of lHarrison et al] 
d2012tl . We take their Lx. boi ,200 and rescale by assum- 
in g a /3-model with Vc es timated using the rc-Lx relation 
of IBohringer et ID d2000l'l and /3=0.67, then correcting to 
Lx.boi ,500 using Eq. |3l The optical luminosities provided are 
calculated for r = 0.5r2oo ~ riooo. By the reasoning de¬ 
scribed previously, this luminosity is corrected to Lr, 500 us¬ 
ing the relation M500 / Miooo oc 1.3. Because the magnitudes 
of the BCG were not reco rded, we rescale all of the optical 
light for these objects. The lHarrison et^ d2012ll are also 
used, and we assign a 0.1 dex error to these values during 
our fit of the fossil scaling relations. 

With the above data sets, we have enough data to as¬ 
semble and quantitatively compare the Lx-Tx, Lx-it„, Lx- 
Lr, Ty-av scaling relations for a sample of fossils and a 
control sample of normal groups and clusters. We do not in¬ 
vestigate the Tx-Lr relation due to the small subsample of 
our control population with both Tx and Lr measurements. 

We fit the equation 


log(y) =a-|-folog(X) 


(5) 


to the data usin g the BCES orthogonal method 
(lAkritas fc Bershadvl Il996ll which accounts for mea¬ 
surement errors in the data as well as intrinsic scatter in 
the fitted relation. We choose to compare the fit of the fossil 
sample to a combined sample of groups and clusters (G-l-C) 
in the same parameter range as the fossil sample. For the 
fossil system data set we exclude NGC 6482 from IkpjotI 
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Table 5. Global properties of the ICM-dominated subsample. 


FGS 


p/rsoo 


kT^p 

[keV| 


^ap 

[z©] 


^X,bol,ap 
|10^'* erg s“^] 


04 

1 

o oi+u.iy 

^•0-L_o.19 

0.401“; 

O.UO_Q 19 

3.5’ (0.71 Mpc) 

5 03+^’^^ 

O.UO_0 19 

1.3+0.1 

14* 

1 

c: 26+0-44 

0-211^8 

7 71 4-0.29 
‘ ' -^-0.29 

4.8’ (1.03 Mpc) 

7 71 +0-29 
' •' -^-0.29 

3.8+0.5 

25 

0.98 

Q 02+0-15 

o.yz_Q 15 

0-28l°;°l 

q Qrj + 0.09 

d.»U_o.o9 

8.5’ (0.92 Mpc) 

3.84l®:°« 

2.4+0.2 

26* 

0.84 

q qq + 0.34 
o.oo_Q 30 

0-191^8 

0 70+0 04 
‘ 0>_o 04 

10.3’ (0.85 Mpc) 

0.821®;°® 

1.9+0.3 

27* 

1 

q qrj+0.33 

O.OU_Q 31 

0.l8l°;« 

q qQ+0.16 

0.00_Q 10 

4.3’ (0.80 Mpc) 

q qO + 0.16 
O.OO_0 15 

1.7+0.3 

30* 

1 


0.30l°;°® 

O.UD o.o6 

6.8’ (0.84 Mpc) 

Q r»^+0.06 

3-06Io.o6 

1.9+0.1 


rsoo 


i'X,bol,r500 
[10^^ erg s“^] 


^500 


Note: Z/x bol the unabsorbed X-ray luminosity in the 0.1-100 keV energy range 
Confirmed fossil system 


and XMMXCS J030659.8+000824.9 from iHarrison et al.l 
1 I 2 OI 2 I) as they are clear outliers. 

We cross-checked the results obtained with the BCES 
method wit h the IDL Astronomy library tool LIN- 
MIX_ERR (|Kellvll200it l. a Bayesian fitting method for lin¬ 
ear regression. The plotted scaling relations and BCES fits 
are shown in Fig. [5] and the best-fitting parameters of the 
relations are recorded in Table[6l Uncertainties on the BCES 
best-fitting parameters are estimated using 10,000 bootstrap 
resamplings. For the LINMIX_ERR fits, the quoted val¬ 
ues are the mean and the standard deviation of the posterior 
distributions for the regression parameters. We investigate 
changing the pivot point of the fits, i.e. rescaling the X and 
Y values in Eq. [5]by a constant, but no difference is found 
in the returned fits. 

We find the BCES fits to the fossil sample are consistent 
within error to the combined groups and clusters fit for each 
scaling relation investigated in this work. The LINMIX fos¬ 
sil and non-fossil fits are for the most part consistent within 
Icr; the i/-intercepts of I/x,boi~Tx and the i/-intercepts and 
slopes of Lx-ctv are consistent within 2cr. These slight dis¬ 
crepancies in the LINMIX fits are most likely due to inho¬ 
mogeneities in the data or the small sample size of both the 
fossil and control populations. 

The global properties involved in these scaling rela¬ 
tions: Lx, Tx, Lopt, cr„, are determined predominantly by 
the shape and depth of the potential well, and are thus well- 
documented proxies for the total mass of the system. Addi¬ 
tional important effects that determine the X-ray properties 
of th e ICM include the entropy structure llDonahue et al.l 
I 2 OO 6 II and non-gravitational heating and cooling processes, 
such as can be caused by AGN or mergers. These factors can 
produce dispersions in scaling relations between X-ray and 
optical mass proxies. Finding no difference in the scaling re¬ 
lations between fossil and non-fossil groups and clusters thus 
indicates fossil systems are of similar mass as non-fossils, 
and on the global scale, the combined effect of mass, ICM 
entropy, and non-gravitational processes that have occurred 
in fossil systems are similar to the combined effect of those 
that have occurred in normal groups and clusters. 


8.2 Comparison with previous studies 

Our result that fossils share the same Lx-Tx relation as non- 
foss il groups and clusters is consistent with previou s stud - 
ies (IKPJOTI : IProctor et al.l[201ll : IHarrison et al.ll2012l . iGldli . 
However, the comparison of optical and X-ray properties of 


fossil and non-fossil systems is a contentious issue in the 
literature. 

The Lx-Lr, Lx-at,, Tx-ctv scaling relation fits of our 
analysis show the relations of fossil systems are consistent 
within error to normal groups a nd clusters. Thi s is in good 
agre e ment with the findings of IHarrison et al. I (l2012tl and 
|Gi3.|Gl3 recorded the first quantitative values of their fit 
to the Lx-Lr relation and found no difference between fossil 
systems (Lx oc and a sample of non-fossil clusters 

(Lx oc While qualitatively we both find no dif¬ 

ference in the Lx-Lr fossil and non-fossil scaling relations, 
there are some numerical differenc es in the returned best¬ 
fitting parameters of our study and lG14l . 

Ou r foss il fit of Lx oc is consistent within 

error to although this is in large part due to the con¬ 
siderable error on both of our slopes. However, our non-fossil 
fit (Lx oc is not within error of the fit determined 

bv lGI4l . Differences in the slopes of our fits could be due to 
mult iple reasons: (1) we use bolometric Lx in our fits, while 
1^ uses Lx.o. i-2.4 kev; (2) our Lx are defined within rsoo 
while the fitted [Gi 3 Lx represent a total luminosity that has 
not been defined within a precise radius; (3) we use different 
fitting methods; (4) we fit our control sample of non-fossils 
over a different parameter space (i.e., one defined to match 
our fossil sample). 

We ch eck to see if we can return more consistent results 
with IG14| by repeating our analysis of the Lx-Lr relation 
using Lx,o.i-2.4keV instead of Lx,boi and expanding the fit 
of our control ‘G+G’ sample to the full parameter space. We 
find the returned fit to the fossil sample (Lx oc 
and to the non-fos sil sa mple (Lx oc are both 

within error of the IG14| fits. And again we emphasize that 
even without the changes m ade h ere, although numerically 
our fits differ from those of [mi. the interpretation is the 
same: fossil systems follow the same Lx-L^ scaling as non¬ 
fossil systems, supporting our conclusion that on the global 
scale, fossil systems have optical and X-ray properties con¬ 
gruent with those of normal groups and clusters. 

Accumulation of multiple differences in data and 
methodology explain the differences i n conclusions between 
our s tudy and those of earlier studies (lKPJ07l : IProctor et al.l 
that find discrepancies in the optical and X-ray scal¬ 
ing relations for fossil and non-fossils. We have compared 
fossil and non-fossil optical luminosities measured from the 
same photometric catalogue and band, avoiding the need to 
make approximative luminosity estimates for comparisons 
between samples. We have also used optical luminosities de¬ 
fined within the same fiducial radius, thus ensuring a more 
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Figur e 5. cry scalin g relatio ns for fossil and no n-fossil s amples. We abbrevi ate this current work as K+. [Zarattini et al.l 

ll2014[l as Z14. | G iraxdi et al.l ll2014f) as Gl4 jMiller et al] ll2012t) as M12. [doctor et al.N201l|) as Pit jKhosroshahi et al.l <l2007h as KPJ07. 
and iHarrison et al.l ||2012|) as H12. The plotted lines are the orthogonal BCES fits to the fossil sample (dashed line) and to the sample 
of groups and clusters (solid line) in the same parameter range as the fossils. 


equal comparison between data pulled from multiple cata¬ 
logues. Additionally, our large sample size of fossils reduces 
the effect of noise to ensure a more reliable comparison be¬ 
tween the fossil and non-fossil samples. 


sets. Furthermore, our results probe the relations of clusters 
and high-mass groups, and consequently it is possible dlf- 
ferences in the scaling relations exist in the low -mass end 
llDesiardins et aDl2014 ; [Khosroshahi et al.ll2014l ). 


We note, however, that our best-fitting parameters for 
both the fossil and non-fossil samples have large errors. 
Thus, a study of fossil scaling relations could be greatly im¬ 
proved in the future by larger and more homogeneous data 
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Table 6. Best fits to the scaling relations. 


Relation (Y-X) 

Sample 

Fitting Procedure 



BCES Orthogonal 

LINMIX 

ERR 



a 

b 

a 

b 

Lx,bol-Tx 

Fossils 

42.48+0.17 

3.21+0.44 

42.49+0.13 

3.39+0.29 


G+G 

42.55+0.09 

3.25+0.14 

42.74+0.07 

3.03+0.11 

Lx.bol-O'u 

Fossils 

30.05+3.60 

4.94+1.29 

28.34+3.22 

5.51+1.14 

G+G 

29.95+1.40 

5.05+0.49 

33.30+0.96 

3.87+0.33 

Lx.bol-Lr 

Fossils 

15.98+3.18 

2.33+0.27 

17.18+2.84 

2.23+0.24 


G+G 

14.47+2.03 

2.45+0.17 

17.05+2.73 

2.24+0.22 

Tx-a-u 

Fossils 

-3.73+2.44 

1.49+0.89 

-4.59+1.67 

1.79+0.59 


G+G 

-3.65+0.44 

1.48+0.15 

-3.92+0.27 

1.58+0.09 


9 SUMMARY AND CONCLUSIONS 

We have presented a detailed study of the X-ray properties 
of 10 candidate fossil galaxy systems using the first pointed 
X-ray observations of these objects. In particular, Suzaku 
XIS data have been used to measure their global X-ray tem¬ 
peratures and luminosities and to estimate the masses of 
these galaxy clusters. We determine 6 of our 10 objects are 
dominated in the X-ray by thermal bremsstrahlung emission 
and thus we are able to measure the global temperatures 
and luminosities of their ICM. This sample of six objects 
has temperatures of 2.8 ^ Tx ^ 5.3 keV, luminosities of 
0.8 X 10"^'* < Lx.boi ^ 7.7 X 10'^'^ erg s“^, and occupies the 
cluster regime in plotted scaling relations. 

Using our newly determined fossil cluster ICM X-ray 
properties, we combine our fossil sample with fossils In the 
literature to construct the largest fossil sample yet assem¬ 
bled. This sample is compared with a literature sample of 
normal groups and clusters, where significant effort has been 
made to homogenize the global Lx, Tx, Lr, and data for 
the fossil and non-fossil samples. 

Plotting the Lx-Tx, Lx-Ut,, Lx-Lr, and Tx-ctv rela¬ 
tions shows no difference between the properties of fossils 
and normal groups and clusters. Furthermore, we provide 
the first fits to three of these relations which reveals the re¬ 
lations of fossils systems agree within error to the relations 
of normal groups and clusters. Our work indicates that on 
the global scale, fossil systems are no different than non¬ 
fossil systems. However, the distinguishing large magnitude 
gap in the bright end of the fossil system luminosity function 
is still unexplained and thus further studies are necessary to 
characterize the properties of these objects and understand 
their nature. 


ACKNOWLEDGEMENTS 

This research has made use of data obtained from the Suzaku 
satellite, a collaborative mission between the space agen¬ 
cies of Japan (JAXA) and the USA (NASA). We thank the 
anonymous referee for valuable comments, K. Hamaguchi 
and K. Pottschmidt at the Suzaku Helpdesk for useful ad¬ 
vice on multiple aspects of our analysis, and D. Eckert for 
helpful discussions and for suggesting the flickering pixels 
issue as an explanation for the excess in the 0.5-0.7 keV 
range. 

Support for this research was provided by NASA Grant 
No. NNX13AE97G, and by the University of Wisconsin- 


Madison Office of the Vice Ghancellor for Research and 
Graduate Education with funding from the Wisconsin 
Alumni Research Foundation. FG acknowledges the finan¬ 
cial contribution from contract PRIN INAF 2012 (‘A unique 
dataset to address the most compelling open questions about 
X-ray galaxy clusters’) and the contract ASI INAF NuS- 
TAR 1/037/12/0. ED gratefully acknowledges the support 
of the Alfred P. Sloan Foundation. MG acknowledges fund¬ 
ing from MIUR PRIN2010-2011 (J91J12000450001). JALA 
has been partly funded from MINEGO AYA2013-43188-P. 
EMG is partially supported by Padua University through 
grants 60A02-4807/12, 60A02-5857/13, 60A02-5833/14, and 
GPDA133894. JMA acknowledges support from the Euro¬ 
pean Research Council Starting Grant (SEDmorph; PI V. 
Wild). 


REFERENCES 

Aguerri, J. A. L., Girardi, M., Boschin, W., et al. 2011, A&A, 
527, A143 

Akritas, M. G., & Bershady, M. A. 1996, ApJ, 470, 706 

Anders, E., & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 
53, 197 

Arnaud, M., Pointecouteau, E., & Pratt, G. W. 2005, A&A, 441, 
893 

Bahcall, N. A., & Kulier, A. 2014, MNRAS, 439, 2505 

Barnes, J. E. 1989, Nature, 338, 123 

Bohringer, H., Voges, W., Huchra, J. P., et al. 2000, ApJS, 129, 
435 

Bohringer, H., Schuecker, P., Guzzo, L., et al. 2004, A&A, 425, 
367 

Buote, D. A., Gastaldello, F., Humphrey, P. J., et al. 2007, ApJ, 
664, 123 

Gavaliere, A., & Fusco-Femiano, R. 1976, A&A, 49, 137 

—. 1978, A&A, 70, 677 

Gravens, T. E. 2000, ApJ, 532, L153 

Gypriano, E. S., Mendes de Oliveira, G. L., & Sodre, Jr., L. 2006, 
AJ, 132, 514 

Dariush, A., Khosroshahi, H. G., Ponman, T. J., et al. 2007, 
MNRAS, 382, 433 

Dariush, A. A., Raychaudhury, S., Ponman, T. J., et al. 2010, 
MNRAS, 405, 1873 

Desjardins, T. D., Gallagher, S. G., Hornschemeier, A. E., et al. 
2014, ApJ, 790, 132 

Doe, S., Nguyen, D., Stawarz, C., et al. 2007, in Astronomical 
Society of the Pacific Gonference Series, Vol. 376, Astronomical 
Data Analysis Software and Systems XVI, ed. R. A. Shaw, 
F. Hill, & D. J. Bell, 543 

Donahue, M., Horner, D. J., Cavagnolo, K. W., & Voit, G. M. 
2006, ApJ, 643, 730 


© 2015 RAS, MNRAS OOP.IRflTl 











Fossil group origins - VI 15 


D’Onghia, E., Sommer-Larsen, J., Romeo, A. D., et al. 2005, 
ApJ, 630, L109 

Eckmiller, H. J., Hudson, D. S., &; Reiprich, T. H. 2011, A&A, 
535, A105 

Eisenstein, D. J., Annis, J., Gunn, J. E., et al. 2001, AJ, 122, 
2267 

Ettori, S., Gastaldello, F., Leccardi, A., et al. 2010, A&A, 524, 
A68 

Fabian, A. C., & Barcons, X. 1992, ARA&A, 30, 429 
Fujimoto, R., Mitsuda, K., Mccammon, D., et al. 2007, PASJ, 
59, 133 

Girardi, M., Giuricin, G., Mardirossian, F., Mezzetti, M., Sz 
Boschin, W. 1998, ApJ, 505, 74 

Girardi, M., Manzato, P., Mezzetti, M., Giuricin, G., &; Limboz, 
F. 2002, ApJ, 569, 720 

Girardi, M., &: Mezzetti, M. 2001, ApJ, 548, 79 
Girardi, M., Aguerri, J. A. L., De Grandi, S., et al. 2014, A&A, 
565, A115 (G14) 

Harrison, C. D., Miller, C. J., Richards, J. W., et al. 2012, ApJ, 
752, 12 

Hess, K. M., Wilcots, E. M., & Hartwick, V. L. 2012, AJ, 144, 
48 

Hogg, D. W. 1999, preprint (astro-ph/9905116) 

Hudson, D. S., Mittal, R., Reiprich, T. H., et al. 2010, A&A, 
513, A37 

Ishibashi, W., & Gourvoisier, T. J.-L. 2010, A&A, 512, A58 
Ishisaki, Y., Maeda, Y., Fujimoto, R., et al. 2007, PASJ, 59, 113 
Jones, L. R., Ponman, T. J., & Forbes, D. A. 2000, MNRAS, 
312, 139 

Jones, L. R., Ponman, T. J., Horton, A., et al. 2003, MNRAS, 
343, 627 

Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 2005, 
A&A, 440, 775 

Kawaharada, M., Okabe, N., Umetsu, K., et al. 2010, ApJ, 714, 
423 

Kelly, B. C. 2007, ApJ, 665, 1489 

Khosroshahi, H. G., Gozaliasl, G., Rasmussen, J., et al. 2014, 
MNRAS, 443, 318 

Khosroshahi, H. G., Maughan, B. J., Ponman, T. J., & Jones, 
L. R. 2006, MNRAS, 369, 1211 

Khosroshahi, H. G., Ponman, T. J., &: Jones, L. R. 2007, MN¬ 
RAS, 377, 595 (KPJ07) 

Kuntz, K. D., & Snowden, S. L. 2000, ApJ, 543, 195 
Kushino, A., Ishisaki, Y., Morita, U., et al. 2002, PASJ, 54, 327 
Lovisari, L., Reiprich, T. H., & Schellenberger, G. 2015, A&A, 
573, A118 

Massaro, E., Giommi, P., Leto, C., et al. 2009, A&A, 495, 691 
Maughan, B. J., Giles, P. A., Randall, S. W., Jones, C., & For¬ 
man, W. R. 2012, MNRAS, 421, 1583 
Mendez-Abreu, J., Aguerri, J. A. L., Barrena, R., et al. 2012, 
A&A, 537, A25 

Miller, E. D., Rykoff, E. S., Dupke, R. A., et al. 2012, ApJ, 747, 
94 

Milosavljevic, M., Miller, C. J., Furlanetto, S. R., & Gooray, A. 
2006, ApJ, 637, L9 

Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 
493 

Osmond, J. P. F., Ponman, T. J. 2004, MNRAS, 350, 1511 
Pointecouteau, E., Arnaud, M., & Pratt, G. W. 2005, A&A, 435, 
1 

Ponman, T. J., Allan, D. J., Jones, L. R., et al. 1994, Nature, 
369, 462 

Ponman, T. J., & Bertram, D. 1993, Nature, 363, 51 
Popesso, P., Biviano, A., Bohringer, H., & Romaniello, M. 2007, 
A&A, 464, 451 

Pratt, G. W., & Arnaud, M. 2005, A&A, 429, 791 
Pratt, G. W., Croston, J. H., Arnaud, M., & Bohringer, H. 2009, 
A&A, 498, 361 


Proctor, R. N., de Oliveira, C. M., Dupke, R., et al. 2011, MN¬ 
RAS, 418, 2054 

Rasmussen, J., &: Ponman, T. J. 2007, MNRAS, 380, 1554 
Santos, W. A., Mendes de Oliveira, C., & Sodre, Jr., L. 2007, 
AJ, 134, 1551 (S07) 

Serlemitsos, P. J., Soong, Y., Ghan, K.-W., et al. 2007, PASJ, 
59, 9 

Stern, J., & Laor, A. 2012, MNRAS, 423, 600 
Sugizaki, M., Kamae, T., & Maeda, Y. 2009, PASJ, 61, 55 
Sun, M., Voit, G. M., Donahue, M., et al. 2009, ApJ, 693, 1142 
Tawa, N., Hayashida, K., Nagai, M., et al. 2008, PASJ, 60, 11 
Uchiyama, Y., Maeda, Y., Ebara, M., et al. 2008, PASJ, 60, 35 
Veron-Cetty, M.-P., & Veron, P. 2010, A&A, 518, AlO 
Vikhlinin, A., Kravtsov, A., Forman, W., et al. 2006, ApJ, 640, 
691 

Voevodkin, A., Borozdin, K., Heitmann, K., et al. 2010, ApJ, 
708, 1376 

von Benda-Beckmann, A. M., D’Onghia, E., Gottlober, S., et al. 
2008, MNRAS, 386, 2345 

Wu, X.-P., Xue, Y.-J., & Fang, L.-Z. 1999, ApJ, 524, 22 
Yoshino, T., Mitsuda, K., Yamasaki, N. Y., et al. 2009, PASJ, 
61, 805 

Zarattini, S., Barrena, R., Girardi, M., et al. 2014, A&A, 565, 
A116 (Z14) 

Zarattini, S., Aguerri, J. A. L., Sanchez-Janssen, R., et al. 2015, 
A&A, 581, A16 

Zhang, Y.-Y., Andernach, H., Caretta, C. A., et al. 2011, A&A, 
526, A105 


APPENDIX A: TESTING FGS24 FOR SWCX 
CONTAMINATION 

The NASA WIND-SWE proton flux light curve displays el¬ 
evated flux levels greater than 4 x 10^ cm“^ s“^ during a 
significant portion of the FGS24 observation (Fig. lAl|l which 
indicates SWCX photons may contaminate the lower E <1 
keV region of the spectrum (see Section To test for 

evidence of this contamination, we repeat the spectral anal¬ 
ysis of Section f6.II for the time intervals where the flux was 
less than 4 x 10® cm“^ s“^. These results are recorded In 
Table IAT] and we find these results are consistent within er¬ 
ror with those of using the full timespan of the observation 
(Table HI). 


APPENDIX B: CHARACTERIZING THE 
SUZAKU XRT PSF 

We determine a radial model for the Suzaku XRT PSF to 
complete our image analysis in Section r5.2l Our PSF charac¬ 
terization employs archival observations of the X-ray point 
source SS Cyg observed for an effective 52 ks between 2005 
November 18 and 19 [Suzaku sequence number 400007010). 
We clean the SS Cyg event files following the same procedure 
applied to our Suzaku observations (see Section |3|. 

The PSF is characterized using the radial profile of the 
stacked XIS0+XIS1+XIS3 image of SS Cyg that has been 
extracted in the 0.5-10 keV energy range and normalized to 
1 (Fig. IBlIl . The average PSF full width at half-maximum 
(FWHM) is found to be ~ 35 arcsec. Our PSF model con- 
sists of the sum of two exponentials, as recommended by 
ISugizaki et all (l2009l l. and thus the model fit to the SS Cyg 
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Table Al. Best fit spectral parameters during a low proton flux time interval for FGS24. 
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Figure Al. Top: the observed XISl light curve for FGS24. Bot¬ 
tom: the WIND-SWE proton flux light curve plotted for the same 
time span. Proton flux has been found to be correlated to SWGX. 
The elevated proton flux levels during the FGS24 observation may 
potentially cause significant SWGX contaminating emission. 


Table Bl. Best-fitting model to the radial brightness profile of 
SS Cyg. 


Component 

Parameter 

Value 

Units 

expl 

Al 

Cl 

^^0,1 

U.^D_q 45 

9 c+0.01 

0 0+115.4 

counts arcsec”^ 

10 ^ arcsec”! 

arcsec 

exp 2 

A 2 

C2 

CO,2 

U.0‘4_q 54 

« 0+38.9 

O.U i74 y 

counts arcsec”^ 

10 ^ arcsec”! 

arcsec 

background 

k 

Q O+’l'l 

10 ^ counts arcsec”^ 


xVd.o.fixi) 405/53 (7.6) 


brightness profile is: 

S{r) = 

where the constant k accounts for the background. The best¬ 
fitting parameters for this model are recorded in Table EH 


APPENDIX C: NOTES ON THE SAMPLE 

FGS03 is a IZ14l verified fossil system. The AGN (2MASX 
J07524421+4556576) associated with the BCG of this sys¬ 
tem i s both confirmed in the optical dVeron-Cettv fc VerorJ 
l2nilt l and radio. The radio emission from this object con¬ 
sists of strong bipolar jets extending 57 arcsec l|Hess et al.l 
l2012|j . This AGN has also been identified as a Type I Seyfert 
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Figure Bl. Stacked and normalized XIS0+XIS1+XIS3 radial 
brightness profile for point-source SS Cyg in the 0.5-10 keV band. 
The best-fitting model, consisting of the sum of two exponen¬ 
tials and a background constant, is plotted in solid blue. Compo¬ 
nents of the model are plotted with dashed lines, and residuals 
are plotted as triangles. Best-fitting parameters for the model are 
recorded in Table IbT] 

dStern &: Laod fioid l. and appears to dominate the X-ray 
emission observed from FGS03. The spectrum of this object 
is better fit by a power-law {Xr = 1-02) than a thermal model 
iXr = 1-17), and no improvement in the fit occurs when a 
thermal component is added to the power-law model. Fur¬ 
thermore, our imaging analysis finds a /3-model poo rly de¬ 
scribes the observed surface brightness profile. \Zl4 find a 
velocity dispe rsion of = 259 km s“^, the smallest dis¬ 
persion of the lSOTl catalogue. Such a low velocity dispersion 
is typically associated with a cool ICM temperature, which 
would explain why there appears to be very little thermal 
emission when compared to a very bright AGN. 

FGSO 4 is a fossil candidate and has the coolest mea¬ 
sured IGM of our sample (Tx = 2.81 keV). The BGG 
of this system contains the blazar NVSS J080730+340042 
dMassaro et al.l l2009l l and in the radio, iHess et al] ll2012l l 
find bipolar jets originating from this source. We do not 
see evidence of contribution from this object in the spectral 
analysis - the spectrum of FGS04 is fit significantly better 
by a thermal model than a power-law (compare a Xr of 1.14 
to 1.43). 

FGS09 is a fossil candidate system at z = 0.125. 
A background z = 0.73 AGN (QSO B1040+0110; 

RA=10:43:03.84, Dec.=+00:54:20.42) is located 15 arcsec 
from the peak X-ray coo rdinates of FGS09. This AG N is 
confirmed in the optical IVeron-Cettv fc Veronll2010ll and 
the radio (iHess et alJ I 2 OI 2 II bands. Based on our surface 
brightness profile and spectral analyses, this AGN is signifi¬ 
cantly contributing to the observed projected X-ray emission 
of FGS09. A large reduced chi-squared of Xr=5-7 is found 
for the /3-model fit to the radial brightness profile. And, a 
power-law model {Xr = 0.92) fits the spectrum of FGS09 
much better than the thermal model {Xr ~ 1.08). 

FGSI 4 is a confirmed fossil system and is the largest. 


hottest, and most X-ray luminous cluster in our sample, 
with rs oo = 1 Mpc, 3x = 5.3 keV, and Lx = 7.7 x 10“*"^ 
erg s“^. iHess et ^ ll2012ll detected radio-loud emission from 
two central sources; however, we did not see evidence of X- 
ray bright non-thermal emission in our spectr al tes ts. 

FGS15 is a rejected fossil candidate llZ14ll . There 
are a number of contaminating sources in the XIS 
FOV of this source. A radio-loud AGN with an asym- 
metric jet is ass ociated with the BCG of this system 
dHess et al.l l2012l l. Within 40 arcsec of the peak sys¬ 
tem X-ray, the background (z = 0.45) quasar [VV2010] 
J1148 03.2+565411 has been identified optically an d in the 
radio (IVeron-Cettv fc Veron|[201(]l : lHess et al.ir2012ll . Of the 
two visually distinguishable point sources excluded In our 
analysis, the object closest to the centre of the system is spa¬ 
tially consistent with the QSO [VV20101 J11 4755.9+564948 
at 2: = 4.32 dVeron-Cettv fc VeronI I 2 OI 0 II . The further 
south removed point source is located at (RA=11:48:08.38, 
Dec.=+56:48:18.64). The closest known spatial match to 
this object is the radio source NVSS J114838+565327 lo¬ 
cated ^2 arcmin away. Our surface brightness profile analy¬ 
sis reveals that a /3-model (+^=5.2) poorly fits the observed 
emission, and additionally the best-fitting spectral model of 
FGS15 is a power-law. For this object, it is possible multiple 
AGN are contr ibuting to the observed emission; however, as 
noted by |Z14| . FGS15 could also be a filament due to its 
small number of constituent galaxies with large differences 
in velocity. 

FGS 24 is a rejected fossil candidate. No associated 
AGN were identified in the literature. However, the spec¬ 
trum of FGS24 is better fit by a power-law than a thermal 
model (compare a Xr of 1.33 to 1.38). FGS24 was observed 
during a period of potentially strong SWCX emission. While 
we found the best-fitting spectral parameters of the full ob¬ 
servation match those of the isolated time interval of low 
proton flux, it is possible SWGX contamination is occur¬ 
ring even during this interval, obscuring the emission from 
FGS24. 

FGS25 is a non-fossil galaxy cluster (Ha). It is the 
second hottest cluster in our sample with Lx = 3.92 keV and 
a corresponding estimated mass of M 500 = 2.4 x 10^^ Mq. 

I Hess et akl (l2012fl find a radio-loud central point source in 
this cluster; however, our spectral analysis indicates no point 
source contribution as the FGS25 spectrum is much better 
described by a thermal model {Xr = 0.96) than a power-law 
model {Xr = 1. 26). 

FGS26 is a lZ14l confirmed fossil with Lx = 3.3 keV and 
Lx = 0.8 X 10“*^ erg s“^. We find no associated significant 
non-thermal signatures in the spectrum. 

FGS27 is a confirmed fossil with measured global prop¬ 
erties of Tx = 3.3 keV and Lx = 3.4 x 10"^"^ erg s“^. Our 
spectral analysis does not indicate contribution of significant 
non-thermal emission. 

FGS30 is a confirmed fossil with measured global prop¬ 
erties of Lx = 3.4 keV and Lx = 3.06 x 10"^“* erg s“^. A radio- 
loud AGN (2MASX J17181198+563 9563) is associated with 
its bright central galaxy dHess et al. |[2 oT|)- The spectrum of 
FGS30 is better described by the thermal model {Xr = 1.05) 
In comparison to the power-law model (Xr = 1.41). 
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